#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include "solve.h"

#include "CGSolver.h"
#include "FDmatrix.h"



void solvefunc(int iter_z,double dt,double dz, double material[][3],double phi[],double k0[],
               double Hpicar[],double H0[],int solvecase,double rhomx[][3]){
    int lda = 3;
    double A[lda*iter_z];
    double b[iter_z];

    int i;

    ///Finite difference matrix
    //PETER@
    //Hier tritt meiner Meinung nach der Fehler auf, ich glaube die Übergabe der Arrays Hpicar und H0funktioniert nicht
    // Muss ich die als Pointer weitergeben?
    createFDmatrix(iter_z,dt,dz,material,phi,k0,Hpicar,H0,rhomx,A,b,lda);
    for(i = 0;i<15;i++){
        printf(" A %f    \n",A[i]);
    }
    /*switch (solvecase){

    case 1:  // CGSolver
        cgmethod(A,b,Hpicar,iter_z,iter_z,lda);
        break;
    case 2: // solve lapack
        //clapack_dgbsv
        break;
    default:
        break;
    }*/

}
